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Abstract 

We present a preliminary measurement of the angle 0i of the CKM Unitarity Triangle using time- 
dependent analysis oi D ^ K^'k~^t:~ decays produced in neutral B meson decay to a neutral D 
meson and a light meson [B — > Dh). The method allows one to directly extract the value of 20i. 
The ambiguity between 20i and tt — in the measurement of sin can then be resolved. We 
obtain = (16 it 21 it 11)°. The 95% CL region including systematic error is —30° < 4>i < 62°, 
thus ruling out the second solution from sin 2^i = 0.726 ± 0.037 at the 97% CL. 

PACS numbers: ll.30.Er, 12.15.Hh, 13.25.Hw, 14.40.Nd 
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FIG. 1: Diagram for the dominant colour-suppressed amplitude for — > Dtt^. 



Precise determinations of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements [l| 
are important to check the consistency of the Standard Model and search for new physics. 
The value of sin20i, where (pi is one of the angles of the Unitarity Triangle [3], is now 
measured with high precision: sin20i = 0.726 ± 0.037 

A new technique based on the analysis of — > DlKgn^ii^jh^ has been suggested re- 
cently Here we use to denote a light neutral meson, such as 7r°, r] or u. The neutral 
D meson is reconstructed in the Kgir^Tc" mode, whose amplitude content is well known. 

Consider a neutral B meson, which is known to be a 5° at time ttag (for experiments 
operating at the T(45') resonance, such knowledge is provided by tagging the flavour of the 



other B meson in the T(4S') 
the B meson is given by 



BB event). At another time, tgig, the amplitude content of 



-\At\/2T 



bO 



S°\cos(AmAt/2) -i- sin(AmAt/2) 

q 



where At = tsig— ^tag, tbo is the average lifetime of the 5° meson. Am, p and q are parameters 
of B^-B^ mixing (Am gives the frequency of B^-B^ oscillations, while the eigenstates of the 



B' 



Here we have 



effective Hamiltonian in the B -B system are \B±) = p\B ) ± q 
assumed CPT invariance and neglected terms related to the lifetime difference of neutral B 
mesons. 

The B ^ Dh decay is dominated by the CKM favored h — > cud diagram as shown in 
Fig. ^ with roughly a 2% contribution from the CKM suppressed h ucd diagram. Ignoring 
the latter, a neutral D meson produced in a B^ decay is a while that produced in a 
B^ decay is a D^. The D meson state produced at time At is then given by the following 
admixture: 



D^oi^t)) = \d°) cos{AmAt/2) - 



D°\sin(AmAt/2), 



(2) 



where we use to denote the CP eigenvalue of h^, and / gives the orbital angular momen- 
tum in the Dh'^ system jsf. 

The next step is to include the matrix element for the decay D K'^tc'^tt'. We follow 
and describe the amplitude for a decay to this final state as f{m\,rn?_), where and 



m 



are the squares of the two-body invariant masses of the K^tt^ 



and -ft'gvr" 



combinations. 



Assuming no CP violation in the neutral D meson system, the amplitude for a D decay is 
then given by f{irn?_,irn?j^). In the Standard Model, \q/p\ = 1 to a good approximation, and. 
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in the usual phase convention, arg(g/p) = 20i. We then obtain, 

MBo(At) = /(m^,m^)cos(AmAt/2) -ie-^2'^i^/,o(-l)7(m^,m^)sin(AmAt/2), (3) 
MBo(At) = /(m^,m^)cos(AmAt/2) -ie+^2'^iao(-l)7(m^,m^)sin(AmAt/2). (4) 

The time- dependent Dahtz plot density, p, is defined by 

e-l^*l/^BO 

{1 + g[A(m^, m^) cos(AmAt) + S'(m^, m^) sin(AmAt)]}, 

|/(m^,m^)P - |/(m^,m^)p 
|/(m?_, m^)p + |/(m^, m?_)p ' 
-2^;,o(-l)'/m(/(m2 , m^)/*(m^, )e+^2'^i) 
|/(m2,m^)|2 + |/(m2_,m2)|2 ' 

where the 6-fiavor charge is g = +1 (—1) when the tagging B meson is a (B^)- Thus 
the phase 20i can be extracted from a time- dependent Dalitz plot fit to B^ and B^ data. 
Note that this formulation assumes that there is no direct CP violation in the B decay 
amplitudes. 

This analysis is based on 386 x 10^ BB events collected with the Belle detector at the 
asymmetric energy e+e~ collider 0]. The Belle detector has been described elsewhere 0. 
We reconstruct the decays B^ D^h^ for = 11^,1] and u and S° D*^hP for = vr^ 
and 7]. We use the subdecays D*^ D°7r°, D° K^n'^Tr', Kg tt+tt", 7r° 77, 
rj — > 77, 7r"'"7r~7r'' and u tt+tt'tt'^. 

Charged tracks are selected with a set of requirements based on the average hit residual 
and impact parameter relative to the interaction point (IP). A transverse momentum of 
at least 0.1 GeV/c is required for each track in order to reduce the combinatorial back- 
ground. All charged tracks that are not positively identified as electrons are treated as pion 
candidates. 

Neutral kaons are reconstructed via the decay Kg tt+tt^ with no PID requirements for 
the daughter pions. The two-pion invariant mass is required to be within 9 MeV/c^ (~ 3a) 
of the K^ mass and the displacement of the tt+tt" vertex from the IP in the transverse 
(r-ip) plane is required to be between 0.2 cm and 20 cm. The direction of a vector in the 
r — if plane from the IP to the tt+tt^ vertex is required to agree within 0.2 radians with the 
combined momentum of the two pions. 

Photon candidates are selected from calorimeter showers not associated with charged 
tracks. An energy deposition of at least 50 MeV and a photon-like shape are required for 
each candidate. A pair of photons with an invariant mass within 12 MeV/c^ {2.5a) of the 
7r° mass is considered as a vr'' candidate. 

We reconstruct neutral D mesons in the i^^vr^vr" decay channel and require the invariant 
mass to be within 15 MeV/c^ (2.5cr) of the nominal mass. D*^ candidates are recon- 
structed in the D^n^ decay channel. The mass difference between D*^ and candidates 
is required to be within 3 MeV/c^ of the expected value (~ 3a). u candidates are re- 
constructed in the vr+Tr^Tr*^ decay channel. Their invariant mass is required to be within 
20 MeV/c^ (2.5 F) of the uj mass. We define the angle 9^ between the normal to the uj decay 
plane and opposite of B direction in the rest frame of uj and require | cos^^;! > 0.3. We 
reconstruct rj candidates in 77 and tt+tt^tt^ final states and require the invariant mass to be 
within 10 and 30 MeV/c^ (2.5cr) of the 77 mass, respectively. The photon energy threshold 



p(m^,m^,At) = 
A(m^,m^) = 
S'(m^,m^) = 
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FIG. 2: From left to right: AE (top) and Mbc (bottom) distributions for the decays to Dir^, 
Duj, Dr] and D*tt^, D*r]. Histo grams represent the data and curves shows the results of the fit. 



for the prompt tt*^ and rj candidates is increased to 200 MeV in order to reduce combinatorial 
background. We remove 1] candidates if either of the daughter photons can be combined 
with any other photon with > 100 MeV to form a vr" candidate. 

We combine and = {tt^, uj, t]} candidates to form B mesons. B candidates are also 
reconstructed from combinations of Z)*° and tt^ or rj candidates. Signal candidates are iden- 
tified by their energy difference in the center of mass (CM) system, AE = {J2i Ei) — -Ebeam, 

and the beam constrained mass, Mbc = \/ E'^ea.m ~ ^iPiY-, where -Ebeam is the beam energy 
and Pi and Ei are the momenta and energies of the decay products of the B meson in the 
CM frame. We select events with Mbc > 5.2 GeV/c^ and \AE\ < 0.3 GeV, and define the 
signal region to be 5.272 GeV/c^ < Mbc < 5.287 GeV/c^, -0.1 GeV < AE < 0.06 GeV (7r°, 
rj — i> 77) \AE\ < 0.03 GeV (cu, 77 -k^-k^-k^). In cases with more than one candidate in an 
event, the one with and masses closest to the nominal values is chosen. 

To suppress the large combinatorial background dominated by the two-jet-like e"'"e^ — > qq 
continuum process, variables that characterize the event topology are used. We require 
I cos^thrl < 0.80, where 6'thr is the angle between the thrust axis of the B candidate and that 
of the rest of the event. This requirement eliminates 77% of the continuum background and 
retains 78% of the signal. We also construct a Fisher discriminant, JF, which is based on the 
production angle of the B candidate, the angle of the B candidate thrust axis with respect 
to the beam axis, and nine parameters that characterize the momentum flow in the event 
relative to the B candidate thrust axis in the CM frame. We impose a requirement on 
that rejects 67% of the remaining continuum background and retains 83% of the signal. 

Figure |21 shows the AE and Mbc distributions for the events in the signal region. For 
each mode, the AE distribution is fitted with a Gaussian for signal and a linear function 
for background. The Gaussian mean value and width are fixed to the values from MC 
simulation of the signal events. The region AE < —0.1 GeV is excluded from the fit to 
avoid contributions from other B decays, such as -B ^ Dh^{'K) where (vr) denotes a possible 
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TABLE I: Number of events in the signal region (A'tot)) detection efficiency, number of signal 
events from the AE fit (iVsig) and signal purity for the B — > D^*'>h? final states. 



Process 


A^tot 


Efficiency (%) 




Purity 




265 


8.7 


157 ±24 


59% 


Doj 


78 


4.1 


67± 10 


86% 


Dr] 


97 


3.9 


58± 13 


60% 


D*7r^,D*r] 


52 




27± 11 


52% 


Sum 


492 




309 ± 31 


63% 



additional pion. For the Mbc distribution fit we use the sum of a signal Gaussian and 
an empirical background function with a kinematic threshold, with a shape fixed from the 
analysis of the off-resonance data. The signal yield is obtained from the fit to the AE 
distribution. The results of these fits to data are summarized in Table HI 

The signal B meson decay vertex is reconstructed using the D trajectory and an IP 
constraint. The tagging B vertex position is obtained with the IP constraint and with 
well reconstructed tracks that are not assigned to the signal B candidate. The algorithm 
is described in detail elsewhere The time difference between signal and tagging B 
candidates are calculated using At = Az/j^c and Az = zcp — ^^tag- 

The proper-time interval resolution function i?sig(At) is formed by convolving four com- 
ponents: the detector resolutions for zcp and Ztag, the shift in the Ztag vertex position due 
to secondary tracks originating from charmed particle decays, and the kinematic approxi- 
mation that the B mesons are at rest in the CM frame A small component of broad 
outliers in the Az distribution, caused by mis-reconstruction, is represented by a Gaussian 
function. 

Charged leptons, pions, kaons, and A baryons that are not associated with a reconstructed 
B^ D[Kg7i^7i^]h^ decay are used to identify the 6-fiavor of the accompanying B meson. 
The tagging algorithm is described in detail elsewhere We use two parameters, q and 
r, to represent the tagging information. The first, q, has the discrete value +1 (—1) when 
the tag-side B meson is more likely to be a (B^). The parameter r corresponds to an 
event-by-event flavor-tagging dilution that ranges from r = for no flavor discrimination to 
r = 1 for an unambiguous flavor assignment. 

We perform an unbinned time-dependent Dalitz plot fit. The negative logarithm of the 
unbinned likelihood function is minimized: 



- 21ogL = -2 



n 

E 



] logp(m^j, m^j, Atj) — log / p{rni^,m^_, At)dm\drn^_dAt 
1 



(6) 



where n is the number of events, m\^, m?_^ and Ati are the measured invariant masses 
of the D daughters, and the time difference between signal and tagging B meson decays, 
respectively. The function p(m^,m^. At) is the time-dependent Dalitz plot density, which 
is calculated according to Eq. © and incorporates reconstruction efficiency, flavor-tagging 
efficiency, wrong tagging probability, background and At resolution. 

We describe the background by the sum of four components: a) B decays containing 
real D mesons, b) B decays with combinatorial D mesons, c) qq events containing real D 
mesons, d) qq events with combinatorial D mesons. The Dalitz plot is described by the 
function f{m'^, w?_) for a) and c) and by the sum of phase space and K*(892) contributions 
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FIG. 3: Dalitz plot distribution for the Dh^ candidates from B signal box (left) and Mbc sideband, 
denote the square of K^tt^ invariant mass. 
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FIG. 4: Kgvr^ (left) and vr+vr (right) invariant mass distributions for the Dh^ candidates. Open 
histograms correspond to the B signal region, hatched ones — to the Mbc sideband. 



for b) and d). The PDF for b) and d) is obtained from an analysis of the events in the 
Mbc-Ai? sidebands. The At distribution for the B decay backgrounds are described by an 
exponential convolved with the detector resolution. For the qq background, a triple Gaussian 
form is used, which is obtained from events with | cos ^thrust I > 0.8. We use the experimental 
data and generic MC to fix the fractions of background components. Figure El shows the 
Dalitz plot distributions for the signal and background events, integrated over the entire 
At range and and B^ combined. Projections of the signal candidate distribution on the 
Mi^Kgii^) (two entries per event) and M(7r^7r~) axes are shown in Fig. EJ 

The procedure for the At measurement is tested by extracting tb+ using 5+ — > 
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FIG. 5: MC invariant mass distributions of D decay daughters produced in the B — > 
D[KgTT^ '7T~]h^ decay chain. Events are generated with 2(/)i = 47°. Only events with good tag- 
ging quahty, r > 0.5, are shown. The dashed histograms show distributions for events with 
q ■ At < —Tgo/2, the sohd histograms show those for events with q ■ At > tqo/2. This range is 
chosen to enhance the visible asymmetry. 



D'^lKgii^n decay. We obtain tb+ = 1.678 ± 0.043 ps (statistical error only), consistent 
with the PDG [ill value 1.638 ± 0.011 ps. 

The potential accuracy of the 0i determination is estimated using a Monte Carlo study. 
We generate D[Kg7i'~^n~]h^ decays using 20i = 47° and process the events with 

detector simulation, reconstruction, flavor tagging and the CP flt. In Fig. El we show the 
invariant mass distributions of the D decay daughters for MC events with q ■ At greater 
than tbo/2 and for events with q ■ At less than — r^o /2. This range is chosen to enhance the 
visible asymmetry. Events with |At| < Tbo/2 are not shown. The MC statistics correspond 
to about 30 times the size of the data. We see clear differences in the two invariant mass 
distributions due to CP violation. Figure El shows the corresponding distributions for the 
data. 

We divide the = [0° : 180°] range into 18 points in steps of 10°. For each point 



9 



cv. 20 
o 

> 18 

0) 

<3 16 



o 
c 

gi 

LU 



14 
12 
10 
8 
6 
4 
2 



_l_l_L 



0.6 0.8 1 1.2 1.4 1.6 1.8 
M(Ksn*), GeV/c^ 



14 
I 12 



o 
« 
c 

gi 

LU 



10 
8 

6 
4 
2 



_L 



_L 



_L 



0.6 0.8 



1.2 1.4 1.6 
M(KsJi"), GeV/c 



1.8 

2 



O 

O 
d 

c 




0.8 1 1.2 1.4 
M(n\'), GeV/c^ 



FIG. 6: Invariant mass distributions of D decay daughters for the B D[Kg'!T'^TT~]hP candi- 
dates. Only events with good tagging quahty, r > 0.5, arc shown. The dashed histograms show 
distributions for events with q ■ At < —tbo/2, the solid histograms show those for events with 
q- At> Tbo/2. 



we perform 30 pseudo-experiments with data samples consisting of 300 reconstructed D7r° 
events. We add 180 background events to each sample, where the background is modelled by 
the shape used in the data analysis. For each pseudo-experiment, we perform an unbinned 
time-dependent Dalitz plot fit. 

Thus, for each input value of 0i we obtain fitted results from 30 pseudo-experiments. 
Prom the means and widths of the distributions of these results we obtain average 0i fit 
results and estimates of their statistical uncertainties. We find the fit results are in good 
agreement with the input values, and the expected uncertainty on 0i is approximately 21°. 

We also study larger ensembles of pseudo-experiments for two 0i input values: 23.5° and 
66.5°, which correspond to sin20i = 0.726 ± 0.037. We verify that the fit value is unbiased. 
The MC pseudo experiments show that the nominal errors from the fit are underestimated 
at the current level of statistics. We obtain a (/>! error of 21° from the RMS width of the 
distribution of fit results, while the errors from the fit have an average value of 15°. With 
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TABLE II: Fit results for the data. The statistical errors shown here are determined from the MC 
pseudo experiments. 



Final state 


01 fit result, ° 


Dtt'> 
Doj, Dr] 

D*-K^,D*7] 


11 ±26 
28 ±32 

25 ±35 


Simultaneous fit 


16± 21 



higher statistics, corresponding to data samples of 1 ab~^ and greater, these two values are 
consistent. Therefore we use the uncertainty from the MC pseudo experiment fit study, 21°, 
as the statistical error for our measurement. The distribution of MC pseudo experiment fit 
results is consistent with, and assumed to be, a Gaussian. This study also shows that, if 
a central value corresponding to one solution of 0i is obtained, the second solution can be 
ruled out at 95% confidence level with our current statistics. 

We have tested for a possible bias in the method due to neglect of the suppressed am- 
plitudes. Due to the smallness of the suppressed amplitude compared to the B^-B^ mixing 
effect, such bias is expected to be small, and indeed we find it to be less than 1%. 

We perform a fit by fixing r^o and Am at the PDG values, using a fixed background 
shape as described above, and using 0i and the background fraction as fitting parameters. 
The de fit results are used to constrain the background fraction. The results are given 
in Table |n] for each of the three final states separately and for the simultaneous fit over 
all modes. Errors are statistical only and determined from the MC pseudo experiments. 
To illustrate the asymmetry, in Fig. [7| we show the raw asymmetry distribution for the 
Dh^ candidates, with an additional constraint to select events consistent with D —>■ Kgp: 
|Mn-+^- — 0.77GeV/c^| < 0.15 GeV/c^. In this case the system behaves as a CP eigenstate, 
with an asymmetry proportional to — sin20i. We do not include D*h^ candidates in Fig. [7| 
since these have the opposite asymmetry. The curve corresponds to the value 0i = 16° 
found in the simultaneous time-dependent Dalitz plot fit over all data. 

The model used for the — > Kg7r~^7T~ decay is one of the main sources of systematic error 
for our analysis. A MC simulation is used to estimate the effects of the model uncertainties. 
We use three models to estimate this error: the first is a simple model from the CLEO 
analysis of KgH^n' decay 3|, the other two are from a similar Belle analysis (gL Esj]. 

Event samples aregenerated according to each model, then we perform a fit with the model 
assumption from [y|. The difference between the input 0i value and fit result does not exceed 
5°. We use different background descriptions to measure the systematic uncertainty due to 
the background parameterization. Several models are used for the Dalitz plot distribution 
of the background: only a uniform distribution and a signal D PDF are compared. For the 
time dependence, we consider cases with only a qq component or only a BB component. 
The difference in 0i results for these models does not exceed 10°. The error due to the 
vertexing and flavor tag is similar to those in other CP analyses [3] (about 1°). 

We have presented a new method to measure the Unitarity Triangle angle 0i using a 
time-dependent amplitude analysis of the D — >• Kg7i~^7i~ decay produced in the processes 
B^ —>■ Dh^. We find 0i = (16 ± 21 ± 11)°. The first error is statistical and determined from 
MC studies. The second is systematic. The 95% CL region including systematic uncertainty 
is —30° < 01 < 62°, thus ruling out the second solution from sin20i = 0.726 ±0.037 at 97% 
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FIG. 7: Raw asymmetry distribution for the D[K^p^]h^ candidates. The smooth curve is the result 
of the fit to the fuh Dahtz plot. 
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